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Domino tilings with diagonal impurities 

Fumihiko Nakano Q and Taizo SadahircJ§ 
Abstract 

This paper studies the dimer model on the dual graph of the square-octagon lattice, which 
can be viewed as the domino tilings with impurities in some sense. In particular, under a 
certain boundary condition, we give an exact formula representing the probability of finding 
an impurity at a given site in a uniformly random dimer configuration in terms of simple 
random walks on the square lattice. 

1 Introduction 

Although the dimer models on planar bipartite lattice graphs have been greatly advanced over 
the last decade (see e.g., [5], [3]), much less is known about non-bipartite cases. This paper deals 
with a non-bipartite lattice T, the dual of the square-octagon lattice. As will be clear later, the 
dimer model on T can be viewed as the domino tiling model containing certain impurities. Our 
main aim in this paper is to study the behavior of these impurities. In particular, under a certain 
boundary condition, we give an exact formula representing the probability of finding an impurity 
at a given site in a uniformly random dimer configuration in terms of the simple random walks 
on the square lattice. 
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Figure 1: Square- Oct agon graph (dashed) and its dual 



We define the dual- square- octagon graph T as follows: the vertices V^(r) is Z 2 which is divided 
into two subsets W = {(x, y) G V(T) \ x + y is even} and B = V(T)\W, and there is an edge 
between v and v' G V(r) if and only if 

v-v'e {±(1,0), ±(o,i)} 
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or 

v,v' G W and v - v' G {±(1, 1), ±(1, -1)}. 

Thus, T is the dual graph of the square- octagon lattice graph (see Figured]). We say a vertex 
is white (resp. black) if it is in W (resp. B). We call an edge connecting two white vertices 
a diagonal edge. The edge set E(T) of Y is divided into two disjoint subsets E\ and E 2 , where 
Ei is the set of diagonal edges and E% = E(T)\Ei. Therefore, the graph F is obtained from the 
ordinary square lattice graph by adding the edges E 2 . We denote by {v,v'} the unoriented edge 
between two vertices v and v'. In the following we sometimes need to orient the edges, and we 
denote by (v,v f ) the oriented edge from v to v'. 

A dimer covering (or perfect matching) M of a graph G = (V(G), E(G)) is a subset of the 
edge set E(G) such that each element of the vertex set V(G) is incident to exactly one element of 
M. We call an edge e in a dimer covering a dimer. We say a subgraph G of T is simply connected, 
if G and T\G are both connected. We say a subgraph of T is normal, if it is simply connected and 
induced by a finite subset of V(r). This paper deals with the dimer coverings of normal subgraphs 
of T. A dimer covering of a normal graph is equivalent to a tilings of the corresponding region by 
square-octagon and octagon-octagon tiles (see Figure [2]). 




Figure 2: Tiling and dimer covering 



For a normal subgraph G of T, we denote 

w G = v(G)nw, b g = v{G) n b. 

Let M be a dimer covering of a normal graph G and let k be the number of diagonal edges in M. 
Then 

2 v ' 

Hence, the number of diagonal edges in a dimer covering M of G is an invariant of G, not depending 
on the choice of M. If a dimer covering M of G does not contain diagonal edge, then it is a dimer 
covering of the ordinary square lattice graph, also known as the domino tiling, which has been 
extensively studied. In this respect, it may be natural to call a dimer e G E 2 of G impurity. Our 
main aim in this paper is to study the behavior of these impurities. In our forthcoming paper 
it will be shown that the local transformations which will be introduced as the t-moves and the 
s-moves in the next section connects all dimer coverings, that is, any dimer covering of a normal 
graph G can be transformed into any other dimer covering of G by applying some sequence of 
the local transformations. This property enables one to construct an ergodic Markov chain whose 
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state space is the dimer coverings. Figure [3] shows the result of a simulation of the Markov chain 
whose stationary distribution is uniform, where we can see that the impurities tend to be located 
near the diagonal edges on the boundary of the graph. 




Figure 3: Markov chain simulation: the initial configuration (left) and the configuration after 10 6 
steps (right) 

The rest of this paper is organized as follows. Section [2] provides some basic properties of the 
dimer model on T. In Section [31 we show an exact result on the easiest case where dimer coverings 
contains exactly one impurity. 

2 Local moves and impurities' orbits 

Let {a, b} and {c, d} be dimers contained in a dimer covering M of a normal graph G, which 
satisfy one of the followings: 

S: a, b, c, d are the four vertices of a unit square. 
T: {a, b}, {b, c} G E 2 and {c, d}, {d, a} e E\. 

In case of S (resp. T), we call the transformation which transforms M into another by replacing 
{{a, b}, {c, d}} with {{b, c}, {d, a}} an s-move (resp. t-move), which is shown in Figure HI 




s-move t-move 
Figure 4: local moves 
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We divide the white vertices W into two parts Wq = TL x 2Z and W\ = Wo + (1, 1) and define 
two graphs A and as follows: A has vertices Wo and it has an edge between v and v' G Wq if and 
only if v — v' G {±(2, 0), ±(0, 2)}, A 1 - is the dual graph of A having vertices W\. Let wi, W2, w^, W4 
be four white vertices which are adjacent to a black vertex b listed in counter-clockwise order as 
shown in Figure El Then one of the two sets {wi,w 3 } and {1^2,^4} is contained in W and the 
other is in W\. Let us assume that Wi,w 2 and b are contained in a normal graph G and w 3 and 
W4 are not necessarily contained in G. For a dimer covering M of G we draw an arc centered at 
W\ (resp. W2) which starts at the middle point of the edge {wi,^} an d ends at a point on the 
edge {wi, b} (resp. {u^, b}) if {u^, b} or {1V4, b} (resp. {w\, b} or {w 3 , b}) is contained in M. Then 
a dimer covering M of G defines curves on the plane composed of these arcs, which we call the 
slit-curves. Figure [6] shows an example of slit-curves. 




Figure 6: Left: Slit-curves generated by a dimer covering. Right: Primary forest (solid) and dual 
forest (dashed) 
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When a t-move r transforms a dimer covering M into another dimer covering M' by replacing 
an impurity e with another impurity e', we simply say r transforms e to e' and simply denote as 

e 1 = r(e). 

Proposition 1. A t-move keeps all slit-curves unchanged. There exists a sequence of t-moves 
which transforms an impurity e to another impurity e' G E 2 if and only if e and e' intersects with 
a common slit-curve. 

Proof. From the definition, a t-move clearly keeps slit-curves unchanged. It is clear that an 
impurity e can be transformed to e' by one t-move if and only if e and e' have a common terminal 
vertex and there exists a slit-curve which intersects with both of e and e'. Now the last statement 
can be easily proved by the induction on the length of the portion of the slit-curve from e to e'. □ 

Corollary 1. If a slit-curve does not terminate on a diagonal edge, it does not intersects with 
impurities. For each slit- curve C, there is at most one impurity intersecting with C . 

Proof. Assume that a slit-curve terminates on an edge {t>i,t>2} G E\ as shown in Figured Then 
the dimer {f2,t>3} must be contained in the dimer covering which generate the slit-curves. For 
the sake of contradiction, assume that an impurity intersects with this slit-curve. The impurity 
can not be transformed to {v 1,^3} in the figure since the vertex V3 must be incident to a dimer 
{^2,^3}, which contradicts Proposition HJ The last statement follows from Proposition [1] and the 
fact that a t-move can transform only one impurity. 




Vi v 2 
Figure 7: A slit-curve terminating on an edge in E\ 

□ 

Proposition 2. 1. A slit-curve does not intersects with another slit-curve. 

2. A slit- curve does not form a loop. 

Proof. From the definition of the slit-curves, [I] is obvious. For the sake of the contradiction, assume 
that a slit-curve C form a loop. Without loss of generality we may assume C is the innermost 
loop. Then the subgraph of G induced by the vertices contained inside C is a tree, since otherwise 
a slit-curve exists inside C and it must form a loop. Since the tree T inside C can not have black 
leaves, i.e., black vertices which is incident to only one edge of T and hence T has odd number of 
vertices, n whites and n — 1 blacks. Thus C must intersect with an impurity. Let us denote by 
Be (resp. Wc) the black (resp. white) vertices which are outside of C and adjacent to vertices 
inside C . Then by the induction on the number of white vertices inside C, we have \Bc\ = \Wc\- 
(See Figure 0) Let b G Be- Then there exists exactly one vertex w G Wc such that {b, w} G M. 
Therefore every element of Wc is incident to a dimer outside of C, which is a contradiction. 
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Figure 8: If a slit-curve C formed a loop, there would be an odd number of vertices inside C and 
no vacant vertex around C. 

□ 

Proposition 3. Remove the edges of G which intersect with slit-curves. Then, each of the con- 
nected components of the resulting graph is a tree, that is, the slit-curves determine a spanning 
forest of G . 

Proof. If a connected component of the graph obtained by removing edges intersecting with slit- 
curves contains a loop, then a slit-curve must form a loop, which contradicts Proposition [2j □ 

Each tree in the forest obtained by removing edges intersecting with slit-curves does not contain 
a path of the form like 

J. 

, i.e., a path bended at a black vertex, hence it can be viewed as a 
tree in A or A^, by removing the edges incident to black leaves, i.e., the black vertices each of 
which is incident to exactly one edge in the tree. We call the set of these trees in A (resp. A -1 ) 
the primary forest (resp. dual forest) obtained from M and denote it by F(M) (resp. F (M)). 

Let M.(G) denote the set of dimer coverings of a normal graph G. We introduce a relation ~ 
on Ai(G) as follows: For two dimer coverings Mi, M 2 G Ai(G), we define 

Mi ~ M 2 

if Mi can be transformed into M 2 by applying a sequence of t-moves. Then the relation ~ is 
clearly a equivalence relation and the set of equivalence classes are denoted by A4(G)/ ~. 

3 Configurations with only one impurity 

In this section, we show an exact enumerative formula of dimer coverings of graphs of special 
shapes, each of which has dimer coverings with only one impurity. 

3.1 Temperley bijection 

Let if be a simply connected subgraph of A. We define the dual graph H 1 - of H in the following 
way. H has the vertices V^if -1 ) consisting of vertices corresponding to the faces of H, more 
specifically, 

V(H^) = {(x, y)eWi = (1, 1) + 2Z x 2Z | (ar, y) is in a face of H} U {/*}, 
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where /* is the vertex taken from W\ so that the graph induced by V{H ) in is simply 
connected in A. Let d* be the number of edges U of A which connect /* and other vetices in 
VlH- 1 ). Then, since H is simply connected, d* is at most 3. H has edges ^(if -1 ) each of which 
corresponds to an edge bounding a face of H. Then, by embedding H and H x simultaneously 
into the plane so that each edge e of H 1 - crosses the corresponding edge e of H only once at the 
middle point of e, where we add a new black vertex, and we obtain a bipartite graph N' . 




H and H ± N G 



Figure 9: By superimposing H and its dual H and removing /* and v*, we obtain a balanced 
bipartite graph N. G is the normal graph induced by V(N) U {f*,v*}. 



Remove from N' the vertex /* and a vertex v* G V(H) which is adjacent to /* in F and 
incident to the outer face of H. Then we obtain a balanced bipartite graph N, which contains 
the same number of black and white vertices. Burton and Pemantle [2] (see also j6]) showed that 
there is a bijection between the set of dimer coverings of N and the set T of spanning trees of H . 
Here, we review this bijection. Let T be a spanning tree of H . Then the edges of H 1 - that do not 
cross the edges of T form a spanning tree of if -1 , called the dual tree and denoted by T x . This 
correspondence makes a bijection between the set of spanning trees of H and that of H . We 
define the root of T (resp. T ) to be /* (resp. v*), and orient T and T 1 - so that they point toward 
the roots. Then the subset M = {{x, \ (x,y) G T or T L } of edges of is a dimer covering 
of N, where (x, y) denotes the oriented edge from x to y. This map T i— > M is the bijection called 
the Temperley bijection [HE]- Conversely, let M be a dimer covering of N. Then the map 



<p: M ^ T = 




{x,^}GM,*G V(H)\. 



is the inverse of this bijection. 

Let G be the normal subgraph of T which is induced by the vertices V(G) = V(N) U {v*, /*}. 
(see Figure [9j) Then a dimer covering M of G contains exactly one impurity. By Corollary [TJ the 
slit-curve which intersects with the impurity terminates at the middle points of the two diagonal 
edges e\ and e\ on the boundary of G, each of which are adjacent to /*. Therefore, there exists 
a unique dimer covering M' such that M ~ M' and e\ G M' . Since M' can be regarded as the 
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dimer covering of N, we obtain the following surjection, 

7T : M(G) 3 M ^ M' G M(N). 
Lemma 1. Let G be the graph described as above and let M be a dimer covering of G. Then 

F(M) = {(poir(M)}. 

LetT 1 - be the dual tree of (f o tt(M) . Then F L {M) can be obtained from T by removing all edges 
incident to f* from T 1 - except for edges in {h, ■ ■ ■ ,ld*}- 

Proof. After removing all edges intersecting with slit-curves from N, the resulting trees whose 
white vertices are in V(A) have no black leaves. Since t-moves keep the slit-curves unchanged, 
the slit-curves can not intersects with T = ip(7r(M)) and T ± , and T is connected, hence F(M) = 
{yjoTr(M)}. □ 

Theorem 1. Let G and H L be graphs as described above. Let Ai be the set of the dimer coverings 
of G and let T be the set of the spanning trees of H . Then the map 



(p : (M/ ~) 9 [M] ^ if o tt(M) G T 



is a bisection. 



Proof. If Mi ~ M 2 then -k{M x ) = tt(M 2 ). Thus ip is well-defined. Since tt is surjective and ip is a 

bijection, (p is also surjective. If <^([Mi]) = tp([M 2 }), then tc(Mi) = 7r(M 2 ) and hence Mi ~ M 2 , 
that is, (p is injective. □ 

3.2 Probability of finding the impurity at a given site 





Figure 10: The impurity intersects with the slit-curve C* which connects the middle points of 
diagonal edges on the boundary. 



Let e = {x,y} G E 2 be a diagonal edge in a normal graph G. Then each class [M] G M.j ~ 
contains at most one dimer covering with impurity e. Thus, to count the dimer coverings with a 
fixed impurity, we can instead count the trees of H corresponding to such dimer coverings. Let 
M G M. be a dimer covering such that T = </?([M]). Then M determines the slit-curves, among 
which the one C* intersecting with the impurity terminates at the middle points of e\ = {f*,v*} 
and e 2 . Thus the slit-curve C* surrounds a tree T* G F (M). (See Figure [101) Therefore, we 
have, 
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Lemma 2. A spanning tree T of H can be represented as T = ipo 7r(M) for some dimer covering 
M containing impurity e = {x,y} if and only if x G T* . 



We choose a spanning tree T of H uniformly at random and define p v by 

p v = Pr(v G T*) (2) 

for each v G V(H). To obtain a uniformly random spanning tree of H, we can instead choose a 
uniformly random spanning tree of H x . By the last half of Lemma [1] and the results of Pemantle 
[3] on the random spanning trees and the loop erased random walks, p v is the probability of a 
simple random walk on H 1 - starting at v to arrive at /* for the first time going through the edges 
in {h, ■ ■ ■ , Id*} ■ By the definition ([2D, it is clear pf* = 1. The probabilities p v 's can be computed 
via the negative Laplacian which is defined as follows. The negative Laplacian A' = (a{j) of H 1 - 
is the ^V r (iJ _L )-dimensional square matrix defined as follows: Each of rows and coluumns of A 1 
corresponds to a vertex of H 1 - and 

— 1 i and j are adjacent, 

% and j are not adjacent, 

4 i = ],i^f\ 

^the number of edges incident to /* % = j — f* '. 

Let A be the matrix obtained by removing the row and column corresponding to /*. Then the 
vector p = ip w )wev\{f*} satisfies 

Ap = b, 

where b = (b w ) we v\{f*} is defined by 



a i,j 



b w 



1 w is adjacent to /* in A J 
otherwise . 



Theorem 2. Let H, H- 1 , A and p v be as described above. Let e = {v,w} be a diagonal edge of 
G m ,n, where v G H^ n . Then the number of dimer coverings containing the impurity e is 

| det A\p v . 

Therefore the total number of dimer coverings of G is 

\M\ = \detA\ 4 Pv + d*-3 \ =\detA\(A(l,A- 1 b) +d* + 1) , 

V veV(R±) J 

where 1 = (1,1,..., 1)' and d* is the number of edges in A 1 - connecting f* and other vertices in 

Proof. By Kirchhoff's Matrix- Tree Theorem (see Chapter 6 of pQ), the number of spanning trees 
of H 1 - is | det A | . Each v ^ f* G Vi^H^) is incident to exactly four diagonal edges of G and /* is 
incident to d* — 1 diagonal edges. □ 
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As an immediate corollary of Theorem [2] we have the following: 



Theorem 3. If we choose a dimer covering of G uniformly at random, the probability of finding 
a given impurity e = {v,w} with v G Vi^H^) is 

Pv 



4£ 



w&V(H ± ) 



Pw + d* 



Example 1. Let H be the graph consisting of n squares in A as shown in the left side of Figure 
[77], and let H 1 - have the vertex f* = (2n + 1,1) and vertices indexed as shown in the middle of 
Figure [771 Then we obtain the graph G as shown in the right side of Figure [771 Then we have 



A 



( 4 -1 
-1 4 -1 
0-14-1 



\ 



4 -1 

V 4 / 

The determinant of A and the probability pj can be evaluated explicitly: 

1 



det A 



2V3 



where A 



± 



2± V3. 



Pj 



Pi 

2V3 



1 - 



A, 



where p x = 2^3 (A+ +1 - A™ +1 ) \ Thus we have 

^ pr_ f \l +1 - A + A™ +1 - A_ 



A_ 



+ 2 



Pi 

2V3 



A + (l + o(l)). 



Thus the probability of finding an impurity e = {j, v} in a random dimer covering is 



1a; ( ^ ) (i + o(i)). 



(0,2) 



(0,0) 



-o- <?- 



(2n,2) 



(2n,0) 



n 
_ 

1 

_ j 



4>-K?>i--^Ti-© 

_ 1 L . 



_ J 




H 



G 



Figure 11: H (left), H 1 - (center), and G (right) 
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Example 2. Let H be the graph as shown in the left side of Figure [7H and let H x be as shown in 
the middle of FigureUM Then we obtain G shown in the right side of Figure\T^ and we have 




A= -1 4 , b=(o,i,i; 



Thus, det A = 56, p = v4 -1 b = (|, |, |)' and the number of dimer coverings is 328. We have 
56 x | = 16 dimer coverings with a fixed impurity incident to the vertex 3, which are shown in 
Figure [73]. 




Figure 13: Dimer coverings with a fixed impurity 
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